www.gusucode.com > 基于matlab编程蒙特卡罗仿真设计源码程序 > 基于matlab编程蒙特卡罗仿真设计源码程序/第四讲 随机数的生成及随机变量抽样/liti7.m
mm=10000 Randnum=unifrnd(0,4,1,mm); xRandnum=zeros(1,mm);yRandnum=zeros(1,mm); for ii=1:mm if Randnum(1,ii)<=1 xRandnum(1,ii)=0; yRandnum(1,ii)=Randnum(1,ii); else if Randnum(1,ii)<=2 xRandnum(1,ii)=Randnum(1,ii)-1; yRandnum(1,ii)=1; else if Randnum(1,ii)<=3 xRandnum(1,ii)=1; yRandnum(1,ii)=1-(Randnum(1,ii)-2); else xRandnum(1,ii)=1-(Randnum(1,ii)-3); yRandnum(1,ii)=0; end end end end xRandnum yRandnum [Y,JJ] =sort(xRandnum) plot(xRandnum(JJ),yRandnum(JJ),'.')